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T/ie motion of a fluid-rigid disc system at 
the zero limit of the rigid disc radius 

Masoumeh Dashti, James C. Robinson 

Abstract 



We consider the two-dimensional motion of the coupled system of a 
viscous incompressible fluid and a rigid disc moving with the fluid, in the 
whole plane. The fluid motion is described by the Navier-Stokes equations 
and the motion of the rigid body by conservation laws of linear and angular 
momentum. We show that, assuming that the rigid disc is not allowed to 
rotate, as the radius of the disc goes to zero, the solution of this system 
' converges, in an appropriate sense, to the solution of the Navier-Stokes 

■ equations describing the motion of only fluid in the whole plane. We also 

prove that the trajectory of the centre of the disc, at the zero limit of its 
radius, coincides with a fluid particle trajectory. 



1. Introduction 



> 

(N 

, 

, The use of rigid tracers for finding the Lagrangian paths of the fluid flow 

■ is based on the assumption that at the zero limit of the rigid body radius, 

fSJ , the trajectory of the rigid tracer coincides with a fluid particle trajectory. 

' Studying the validity of this assumption has been the motivation for the 

problem that we address in this paper. We consider the system of one rigid 
disc moving with the fluid flow in R^, assume the disc does not rotate around 
its centre, and show the convergence of the trajectory of the centre of the 
disc to a fluid particle trajectory. 

We consider the two-dimensional domain occupied by an incompress- 
ible fluid of density 1 and viscosity ly and a rigid disc of radius r and density 
p, which here we set it equal to the density of the fluid. We assume that 
the motion of the fluid is modelled by the Navier-Stokes equations with 
no-slip condition on the boundary of the rigid disc and consider the rigid 
body motion to be described by the equations of the balance of linear and 
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angular momentum: 



— - i^AU + {U ■ V)U + VP = F, 

V • [/ = 0, 
U = h{t)+uj{t){x-h{t))^, 

Mh{t) = - [ Sndr + M^, 



xeR'^\ B{t) (1) 

x e R2 \ Bit) (2) 
X e dBit) (3) 



JdB{t) 



(4) 



Juj{t) 



f (x - h{t))^ ■ Sndr + Jr], 



(5) 



JdB(t) 



MO) -(0,0), 



/i(0) = /lo e K^ w(0) = Wo e M. 



;7(x,0) = [/o(x), 



xeR2\5(0), (6) 
(7) 



In the above system t € [0, T] , i3(t) is the region occupied by the disc at time 
t, U is the velocity vector of the fluid, P the pressure scalar field, h{t) the 
position of the centre of the rigid disc at time i, lo its angular velocity and 
M = pirr^ and J = are its mass and moment of inertia respectively. 

The stress tensor S is defined as 



The body force applied to the fluid is denoted by F, and and Jr] are the 
force and torque applied to the disc. The solvability of the above system, and 
also the equivalent system in a three-dimensional domain, is known. When 
the domain of the motion is the whole of M^, this is shown by Judakov 
(1974) and Scrrc (1987) using Galcrkin approximations. In the case of a 
bounded domain, the global existence of at least one weak solution is proved 
by Hofi'mann & Starovoitov (1999) for the two-dimensional case and by 
Conca, San Martin & Tucsnak (2000) and Gunzburger, Lee & Seregin (2000) 
for three-dimensional domains under some constraints on collision of the 
rigid body and the boundary of the domain. The method of Hoffmann & 
Starovoitov (1999) is based on approximating the rigid body by a very 
viscous fluid, Conca et al (2000) write the system in a coordinate system 
attached to the body and use a Galerkin method, and Gunzburger et al 
(2000) approximate the system by time-discretized problems. 

In their 2004 paper, Takahashi and Tucsnak show that the system of a 
fluid-rigid disc in the whole plane has a unique global strong solution using 
a contraction mapping approach. This result is shown to be true for the case 
of a bounded two-dimensional domain as well by Takahashi (2003). When 
the motion is in M^, the existence of a local unique strong solution is proved 
by Galdi and Silvestre (2002). 

In the case of the motion of several rigid bodies in a fluid, at least one 
weak solution is known to exist (Desjardins & Esteban 1999; Desjardins 



S = -Pld + 2uD{U) 



where 
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& Esteban 2000; Feireisl 2002; Grandmont & Maday 2000; San Martin, 
Starovoitov & Tucsnak 2002). The uniqueness however is not known. 

We also note the papers by Robinson (2004) and Iftimie, Lopes Filho 
& Nussenzveig Lopes (2006). In Robinson (2004) a simphfied model of the 
motion of a fluid-particle system at the limit of zero radius of particles 
is studied. The simplification is in considering only momentum exchange 
between fluid and particle and assuming that the domain of the fluid does 
not change in time; and it is shown that in the limit the motion of the fluid 
is described by the standard Navier-Stokes equations. Iftimie et al (2006) 
consider the two-dimensional asymptotic motion of a fluid outside a small 
obstacle as the size of the obstacle tends to zero. They prove that the limit 
flow satisfies the Navier-Stokes equations in the full plane. 

Knowing that a unique strong solution of ((I])-© exists for any r > 
by the result of Takahashi and Tucsnak (2004), in this paper we study the 
limiting behaviour of the solution of (HJ-© as r — )■ 0. To do this, we do 
not take into account the torque exerted by the fluid on the particle, i.e. 
we do not allow the particle to rotate. Of course, one could incorporate this 
into a model with angular momentum by considering ujq = and imposing 
a torque 



so that a; = 0, but this simply has the effect of removing equation ([5]) 
from all our calculations. It is of course desirable to study the more general 
system when the disc can rotate, but with a nonzero uj we have not been 
able to obtain the required uniform estimates (see Remark [1] in Section 
In order to prove the convergence of the velocity field as r — >■ one requires 
uniform bounds on the fluid velocity field for r > that are not provided 
by the method of Takahashi and Tucsnak. Here, we obtain such uniform 
bounds over a small enough time interval [0,T]. We then show that at the 
zero limit of r, the velocity field satisfies the Navier-Stokes equations in 
the whole plane. Then using the properties of the strong solutions of the 
Navier-Stokes equations, we can prove this convergence for all times. 

In what follows, in Section[2]we recall the result of Takahashi & Tucsnak 
(2004) on the existence of strong solutions, and also their functional setting 
for the problem, which we will follow in this paper. In Section [3] we prove 
some of the properties of the functional spaces introduced in Section [21 
In Section |4] we show that the solution of the above fluid-rigid system is 
bounded independent of the radius of the rigid disc, in the time interval 
[0,T] for a small enough T which is also independent of r. In Section [5] we 
show that as r — > 0, the above fluid rigid system converges to the Navier- 
Stokes equations in the whole of for all times and the trajectory of the 
centre of the disc converges to a fluid particle trajectory. 
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2. Preliminaries 

In this section, following Serre (1987), Conca et al (2000) and Takahashi 
and Tucsnak (2004), we write (HJ-© with respect to a coordinate system 
moving with the disc centre which results in the equations of motion of the 
fluid and rigid body in domains which are fixed with respect to time. We 
will also state and consider the functional setting introduced in Takahashi 
and Tucsnak (2004) for this problem. 

Using the following change of variables 

y{t)=x{t)~h{t), 
u{y, t) = U{y + h{t),t), p{y, t) = P{y + h{t), t), 
a{y,t) = -p{y,t)ld + 2vD{u){y,t) 

the equations of motion can be written as 
du 

— - uAu + [u ■ V)m - (hit) ■ V)w + Vp = 0, 
ot 

V - u = 0, 

u = h{t) +uj{t)y^, 

Mh{t) = - I andr, 

JdB,. 

JCj{t) = - y-^ ■andr + Jrj, 

JdBr- 

u{y,0) = uo{y), 

h{0) = (0, 0), h{0) = /lo e uj{0) = Wo e R. 
with Br = B{0) and 12^ = \ Br- 

We define 

Hr = {ue[L^{R^)]^ -.V ■u = inM^ D{u) = { 
Vr ^ {u e [H^m^)]^ : V ■ u = inR^ D{u) = 

and consider the inner product on Hr to be 

(m, v) — {u ■ v)dy + p {u ■ v)dy 
J Qr Jsr- 

which is equivalent to the standard inner product on [L2(R2)]2. We denote 
the orthogonal projection of [L/^^M?)]"^ onto Hr by P^. 
We let 



y Cz fir 


(8) 


y Cz f^r 


(9) 


y e dilr 


(10) 




(11) 




(12) 


y g Hr, 


(13) 




(14) 


in Br} 


(15) 


in Br} 


(16) 



D 



(Ar) = {u£ [ffi(R2)]2 : u e [H\nr)]^, 

V-M = inR^ D{u) = in (17) 
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and for any u G D{Ar) define 




—vAu in Qr 

^ D{u)nAv + /^^^ ■ D{u)ndy\v^ in 



and 

ArU = frArU. (18) 

It can be shown (Takahashi & Tucsnak 2004) that Ar is self adjoint: 

(ArU, u)l2(r2) ^2v I D{u) : D{v) dy = {u, Arv)L2(jg2y 
J n 

We can extend u e i^(i?r) to a function in iJ^, defined on the whole of R^, 
by letting 

u(y) — h{t) + wj/^, in B^. 

We define 

/3(«-/.,^) = |(("-'^^-^)"^^^'- 
[ in Br. 

Then using Ar and defined as above, equations (l8))- (fT4)) can be written 
as 

— +ArU = -¥rP{u-h,u). (19) 

Existence of a strong solution of ([5])- (|14l) is shown by Takahashi & Tuc- 
snak (2004). We recall their result: 

Theorem 1. (Takahashi & Tucsnak 2004) i) Suppose uq E H^{flr) with 
\/ ■ uq = in ilr and uo{y) = Iiq + ^oj/^ on dfir- Then, there exists a 
unique strong solution of ^- {14% satisfying 

u e L2(0, T; H^{Qr)) n C(0, T; H^iHr)) n ff^(0, T; L'^iflr)), 
p e L^{0,T;H^), he H^iO^T;^^), lu e H\0,T;R), 

for any T > 0. 

ii) Let Uq G Lp'{flr) with V • mq = 0. Then the system ([3j- |J^[ ) /las a 
unique weak solution with u G L^{{),T;H^{Qr))^C{0,T-L^{nr))- 

The estimates in the result of Takahashi and Tucsnak (2004) depend on 
r via ||/i||loc(o,T;R2) (they use M|/ip < c(l + ||Mo||ifi)) and the constant c 
in the inequality || V^M||i2(^2) ^ c||^rWr||L2 used in their proof. But here, 
to study the solution when r — 0, we need uniform bounds on the velocity 
field. In Section 2] we show that if Lo{t) = for < > 0, one can find such 
uniform bounds on a time interval [0,T], with T only depending on the 
initial velocity field uq. But first we prove some of the properties of the 
orthogonal projection which we will need later. 
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3. Properties of the orthogonal projection 

In this section we prove some results characterizing P^, which we will 
need in the following sections. In the first lemma we identify the orthogonal 
complement of H^. Then we use it to show that the image of a sufficiently 
regular and divergence- free function under P^ converges in [L^(R^)]^ to the 
function itself as r — >■ 0, which will be key in proving the convergence of 
the fluid-rigid body system ([TS]) to the Navier-Stokes equations in the whole 
plane in Section [5l We note that the results of this section are in the general 
case where uj{t) can be nonzero. 

Lemma 1. Let 

Hr ^ {u£[L^{«?)f -.V ■u = Q mR^ D{u) = Q in Br}, 

G^ = {ue [L\M^)f : u = Vgi, Vgi G [L\R^)]\ q, e iL(R')}, 

G2^{ue [L2(M2)]2 : V • u = m 

u = V(j2 in with \/q2 G [L^(i7T.)]^, 92 G -^z'ocl'^i-), 

u — 4> in Br with (j) G [L^(i?r)]^, and / ip ■ dy = 0}. 

Then Hr, Gi and G2 are mutually orthogonal and 

[i2(R2)]2 ^ Hr®Gi®G2. 

Proof. The proof is mostly based on the ideas in the proof of Lemma 4 
of Conca et al (2000). Since Hr and G2 are subsets of the set of square 
integrable divergence-free functions, Gi is perpendicular to both of them 
(noting that since we have assumed p = I, the inner product over the above 
spaces is the standard L^-inner product). To show that Hr-LG2, we consider 

V(72 in [2r ^ „ 
<p m Br 

and V £ Hr with v{y) — Vy + oo^y-^ for y £ Br where K, G and G M 
are constant. We have 

{v,w)^ / v-\7q2dy+ / {Vy + uj^y^) ■ <j) dy 

= K- / q2ndy+V,-[ (l)dy. (20) 

JdOr JBr 

Since V?/ is the identity matrix and (p is divergence-free, we can write 
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Now since w G [i^(]R^)]^ is divergence-free, the normal components of 
and V(72 on the boundary, by Theorem 1.2 of Temam 1977, satisfy 



/ {(j)-n-Vq2-n)g{y)Ay 

Jdn-r 



(21) 



for any g € H^^'^{df2r)- Therefore, with g{y) = y and since Aq2 = 0, we 
have 



/ ■ {~Tn)ydy = - / Vq2 • n 
Jdn-r Jan,. 



V • (V(Z2y2] 



dy 



Vy V(72 dy 
/ 9^ 92 % dy = - / 5j (g2 ^ij ) dy 



92 % rij dy = - / (72 ndy 



where i,j = 1,2, 9j = d/dyj and 



1, if i = j, 
0, otherwise. 



Therefore 

/ ^dy^~ [ (Z2ndy (22) 

J Br Jdnr 

and by (EOl), (w, w) = 0. 

Now we consider an arbitrary u G [L^(R^)]^ and show that there exist 
a unique v £ Hr, with 

w(y) K + Wi,y^ for ye Br (23) 

where K € and w^, S R are constants, and also unique Vqi £ Gi and 
w € G2 with w = Vq2 in ■'^r and w — (f> in Br, such that 



u = V + Vqi + w = 
For qi we solve 



u + Vgi + Vq2 in ilr, 

Vr + UJ^y-^ + Vqi + (j) in Br- 



rZ\(Zi(y) = V-u(y), y e R^, 
[ Vqi — > 0, as \y\ — > 00. 



(24) 



The above system has a solution qi E Lj^^(R.'^) with Vgi £ i^(R^), which 
is unique up to a constant. Clearly, Vgi S Gi. 
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We set cj): [L^{Br)Y to be 

4) ^ u-Vqi-Vu - i^JvV^ , (25) 

with and w^, as in (1231) . Integrating the above equation over and using 
(HH) we can find Vy in terms of q2 : 



^- = -^{1 (^i-V<zi)dy+ / gandyV (26) 



Now, for 52 "we write 

Z\g2 = 0, in i?., 



(27) 

Vg2 ■ n = (u — Vgi) • n — 14 • n, on dQr- 
Let 

^''' = {'7eiL(^.):VgG[L2(r2.)F}, 
with the scniinorm ||g||£)i.2 = || Vq||i2(j-2^). For any ^ e D^^^, we have 

= (Z\g2,e) - (Vg2,V0+ / Uu~Vqi)-ndy- f ^V^ndy 
= (Vg2,V0+ / C(«- Vgi) -ndj/ 
\ \ ( 92ndy ) ( / ^ndy 

and therefore 

= [ f(M- Vgi) -ndy- ^ f / {u-Vqi)dy]([ ^ndy 

(28) 

Consider [q] ~ {p e -D^'^ : p ^ q + k for some A: G K} and let D^ '^ be 
the the space of all equivalence classes [q] with q E D^'^. Then D^'^ is a 
Hilbert space (Galdi 1994, Theorem II. 5.1). We now use this fact and the 
Lax-Milgram theorem to show the existence of a unique (up to a constant) 
solution to . Let R = 1 + r and Bji an open ball of radius R and centre 
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(0, 0). For any q £ D ' , we have 



/ 

J a fir- 



qndy 



< \df2, 



|l/2 



< Cir) \\q\\m(f2,.nBn) 

< C(r) (||-7||L2(,,^nSH) + I|V'7|1l2(^^)) 

< C(r) ( ||V<7|U2(^,^) + / qdx +\ 

\ J Br 

< C{r) {\\^q\\LHnr) + \\q\\LHBR)) 



Vg|| 



by the Sobolev trace theorem and the Poincare inequahty. Therefore the 
left-hand side of (^5]) is hnear and bounded. It is also coercive since setting 
12 = S, gives an expression bounded below by || Vi^|p. The right-hand side of 
(PS)) is a linear functional of £, as well. It is also bounded because u — Vgi 
is divergence-free and we can write 

/ C(M-Vgi)-ndy < / • (m - Vgi ) dy 

JdOr Jn,. 

< Vqi||i2(R2) ||VC||l2(j-2^). 

Therefore by the Lax-Milgram Theorem there exists a unique q2 G D^'^ 
(and thus a unique up to a constant q-z G D^-^) that solves 
Having qi and q2, we can find vln^ in terms of u by seting 

V ^ u — Vgi — V(72 in ^r- 

It remains to uniquely define Wi, (of (|23l) ) and (j) in terms of u. We let 

= 4 / ^ • 2/-" dy (29) 

(noting that Jg Vgi • y-*- dy = 0), and this then gives a unique by ((25|) . 

We need to verify that w £ G2 and v € H^- We have V • = in 
and Z\(/2 = in f2r- Therefore, for any e C^(]R^) since Vg2 • n = (/> • n, 
we can write 



wVV'dy^ / u;V7/'dy+ / luVf/'dy 
R2 Jr2^ Jb,- 



Vq2 ■ nipdy — / • n -0 dy = 



implying that V • w = in the sense of distributions. 

In showing v e Hj. also, the only condition which is not obvious is 
V • ti = which, having Vq2 •n = M- n — V^-n — Vgi • n, can be verified in 
a similar way to the proof of V • w = 0. □ 



10 



Masoumeh Dashti, James C. Robinson 



Before proceeding to show an appropriate bound on the i^-norm of the 
difference of a divergence- free iJ^-function and its image under P^, we need 
to prove the following lemma which shows how the constant of the trace 
inequality for the domain i?^ depends on r. We note that in this paper we 
use the result of this lemma only for f2r, whose boundary is a circle, but 
the more general statement seems interesting: 

Lemma 2. Let £r he an exterior domain with a hounded houndary d£r of 
diameter 2r and having the uniform -regularity property. 

Assume that the radius of curvature of d£r at all points is bounded below 
by c\r, with c\ a constant independent ofr, and any x €:'S?\£r is contained 

in an open hall of radius C2r with C2 < ci lying inside \ Sp- 

Then for any positive a < 1 there exists a constant c depending only on 
a such that for all u G H^{£r) 

Proof. Since d£r is compact, a finite open cover {t/j}j^i of d£r exists. We 
consider each Uj to be a ball of radius C2r/2 with its centre on d£r and 
note that since the diameter of £r and the radius of the curvature of d£r 
are appropriately bounded, N is independent of r. 

As in Section 5.21 and 5.22 of Adams (1975) we define 9j to be the 
smooth transformation mapping B = {y & ■.\y\ < 1} onto [/, such that 
Uj n d£r = \!/j{B[)) with Bq = {y E B : 7/2 = 0}. Then for Uj as above and 
letting X = ^{y), since the radius of curvature of d£r is bounded below by 
cir, we have 

dx ' 

I-— ^1 < cr for i,k= 1,2 
dyk 

and 

I d{yi,y2) I ^ c_ 
d{xi,X2) ~ r^' 

with c depending on ci and C2. For each Uj (1 £r, j = 1, ■ ■ ■ , N , there exists 
an extension operator Ej such that for any u G H^{£r), 

Ej [u] = u, in Uj n fr- 
aud 

ll-E^jMllHi(R2) < Kj\\u\\Hi(er) 
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with Kj independent of r (Adams 1975, Theorem 4.32). Therefore we can 
write 

/ Iwpdx < V / \E^[u]{x)\'^ dx 

<crJ2 f |i?,M(<Z',(2/))pd2/ 
] "'-1 

J 

j 
j 

< J2 (cr-" \mu]\\miR-^)+<''r\\DE,[u]\\l,^^,^) 
j 

<cr"||u||2^i(£^) 
noting that N is independent of r. □ 

In what foUows we will also use the Sobolev inequality 

ll^llL'!(r2,) < c\\u\\yy2/p-2/q,p(^Q^^) 

where the constant c depends on p, q and the cone C determining the cone 
property of J?^ (Adams 1975). Since the cone C for the exterior domain J?^ 
does not depend on r, the constant c in the above inequality is independent 
of r. 

In the last lemma of this section we show that if the vector field u is 
divergence-free and regular enough, its projection under is arbitrarily 
close to u in L^(R^), if r is sufhciently small. The result of this lemma is 
essential in showing the convergence of equations (l5))- dTil) to Navier-Stokes 
equations. 

Lemma 3. Let u e [_ff"'^(]R^)]^ be divergence-free and assume that r < 1. 
Then for any positive a < 1 there exists a constant c depending on a such 
that 

||-u - Pr-u||i2(R2) < cr"/^ ||m||^i(r2). 
Proof. By Lemma [U since u is divergence- free Vgi — and we have 



u = VrU + V(72 in 
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where (72 satisfies (P7|. Letting ^ = 92 in we obtain 
l|V<Z2||i2(,,^) + ^ <72ndj/^ 

g2U-ndy ( / udy]( q2ndy 

dn^ T^f^ \Jb,. J \Jdn,^ 

< [ q2U-ndij + [ [ udy] + ( / q2ndy 

and therefore 

2 



1 ' " 



l|Vg2|L.(^.) + ^^^^^<72ndy^ 

</ q2u-ndy+-^U udy\ . (30) 
Dropping the second term in the left-hand side, we can write 

\\^Q2\\Un^) < J^^ S/q2udy + ^ (\\u\\L2^B,.)\Br\'^^y 
< l|Vg2||L2(^^^) ||u||l2(r2) + ^||u||i2(B^), 

which gives 

llVg2||L2(^-2^) < c||u|li2(R2). 

Let 

Or = {xe : dist(x,(9/?r) < 2}. 

If q2 solves (P7)) . 

is a solution as well and therefore we can assume that the average of 92 
is zero in Or- By continuity there exists a line passing through (0,0) and 
dividing Or to two bounded and simply connected domains O}. and such 
that the average of q2 is zero on both of them. Noting that the boundaries 
of Or and have strong local Lipschitz property, there exist two extension 
operators Ei and E2 such that for any u e H^{Or), and i = 1, 2 we have 

E,[u]^u, in 01, 

and 

||£^iMIIhi(r2) < ||u||hi(c).), 

with independent of r (Adams 1975, Theorem 4.32). Therefore by The- 
orem 4.4.4 of Ziemer (1989) we have 

||92|U2(o.) < C(|0;|)||Vg2|U2(o.), ^ = l,2, 
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which imphes that 

lk2|lL2(o^) < c||Vg||L2(f2,) 
< c||M||i2(K2). 

Having this, we go back to ((30)) and write 

2 



u Ay 



< \\l2\\L'^(an^)\W\\L^(dn^) + ^;^II^IIlv(i-»)(b,) I-^'-I 
Arguing similar to the proof of Lemma [2j we can show that 
h2\\L^{dn^) < cr"^^ |k2|Ui(o,) 

and therefore 

l|V92||i2(^,) + ^(/^^^^92ndy) 

< cr"/2 {\\q2\\mo,.) + l|Vg2|lL2(^,^)) \\u\\Hiin^) + cr^" \\u\\ 

< cr"^^ ||u||l2(r2) ||m||//i(r2) + cr^" II"II^i(r2)- 
Hence, we obtain 

||u-P,u||i2(,2^) = ||Vg2||L2(f2.) < cW2||u||^i(R2), 

and 

1 

r 



l+a 



2 



g2ndy 



(31) 



To show the bound on ||u — Vru\\]^2(^g^j, we note that 

Pru{y) = + (^vy^ 

for any y e Br, with K, e and G M defined as in and (1^ . 
Therefore 

Pru||i2(s^) = \\u- (K + WDy^)||L2(i3,) 

< c||M||i2/(i-c)(B^)|B,r/2 + cr\V,\+cr^\uj^l 
using Holder inequality. From ((26|) and (j3ip we have 



r|K|<-||^i||LVa-.)|i?r|^^+"'/' + 



g2nd?/ 



Of?,- 



< cr"||w||//i(R2) + cr"/2 ||w||hi(r2). 
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By we can write 
Therefore 

\\u-FrU\\L2(^B,^) < Cr"/2 ||u||/fl(R2), 

and the result fohows. □ 



4. Local uniform bounds on the fluid velocity field 

Here under the assumption of zero angular velocity for the disc for t > 
0, we obtain uniform in r bounds on the velocity field. To impose this 
assumption in the model (I5))- P^ we set 

wo = and Ji] = {x - h(t))^ ■ SndR (32) 

JdB(t) 

We show that for some small enough time interval after the initial time, 
one can bound the solution of ((8))- (ll4l) uniformly in r, in some appropriate 
spaces. We then prove that these uniform bounds result in strong conver- 
gence of the solution in an appropriate sense. 

To show the uniform bounds, we need the following estimate 

\\^'^Ur\\L^(a^) < c\\ArUr\\L2{f2,^) 

with c independent of r. We first prove this estimate. 
Lemma 4. Let u he the solution of 

Au + Wp^ /, 171 ]R2 \ Br 

V • u = 0, m R2 \ Br (33) 
u = hj for X e dBr 

with f S L^{f2r) and the constant h £ M.'^ . Then, there exists c > inde- 
pendent of r and h such that 

\\^Mmf2^) <c\\f\\mn^). (34) 

Proof. Let y — x/r and Vj^ and Ay the gradient and Laplacian operators 
in y coordinates. Then 

AyU + rS/yP = r^/, in \ Bi 

\'y-u = 0, in ]R2 \ Bi (35) 

u — h, for y e dBi 
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Define 'P{x) — h for any x gM.'^ and let v = u — <P. The vector v satisfies 
AyV + r\/yp^r'^f, inK2\Bi 

Vy • w = 0, in \ Bi (36) 

V ~0, for y e dBi 

Now we use an estimate of Galdi (1994, equation V.4.14): 

\\Dlv\\<cr'\\f\\L2 (37) 

with c a constant independent of r. The above estimate holds if S, the set of 
solutions of when / = contains only u = 0. In our case u € L^{f2) and 
^(y) h as y oo. Therefore u — — /i = o(log |y|) and by Theorem V.3.5 
of Galdi (1994) S = {0}. Hence (|37)) is valid. We now revert the coordinates 
to X in (1371) and the result follows. □ 



Remark 1. It is for the result of the above lemma that we need to assume 
that the rotation of the disc is zero. For a nonzero uj where u = h + xuj/\x\ 
on dBr, we have not been able to obtain a uniform bound on ||V^M||i2(^^). 

Now we can show the uniform estimates: 



Lemma 5. Let i32\) hold and consider {ur,hr) to be a strong solution of 
^-(1^ with Ur{0) e H^{Qr) and hr{0) G M^. Then, there exist T and 
C depending on ||itj.(0)||, ||u,.(0)||/fi and \hr{0)\ but independent of r, such 
that Ur is bounded by C in 

Proof. Taking the inner product of 

— h ArUr + Pr l3(Ur — /ir, Ur) ~ 

at 

with Ur € Hr, we obtain 
1 d 



Since 



2 ' i^||VMr||i2(f-2,_) = - / {{Ur - hr,Ur) ■V)Ur ■ Urdy. 



{{Ur ~ hr) ■ V)u • U — -{{Ur — hr) ■ V)|up = -V • (|ltp(Wr — hr)) , 

have 

{{Ur — h, Ur) ■ \/)Ur ■ Ur dy = — j |WrP(Wr ■ TO. ~ hr ■ u) dy ~ 0. 

Therefore 
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which after intrgration gives the following estimates 

\\Urit)\\L2 <\\Ur{0)\\L2 and / \\Ur\\mm^)dt<c\\Ur{0)\\l2 (38) 

Jo 

for almost every t E [0,r]. We also note that since 

hr ■ n = Ur{x) ■ n for almost every x £ dfir 

we have 

\hr\ < c\\Ur\\w^'i(i2r.)- (39) 

with c independent of r. To show this, for any x e df2r, we consider the 
line Lx tangent to dSlr at x and dividing to two parts. Using the trace 
theorem for as the boundary of the part that does not contain B^, we 
can write 

\Ur{x)\ < c||Wr||ioo(i^) < c\\Ur\\w^.i{nr) 

which implies the previous inequality. 

Now we take the inner product of with ArUr G [i^(R^)]^ and obtain 

< II \Ur\ |VWr| \ArUr\ Wl^Q^) + \K\ \\ |Vlir| \ArUr\ \\L^(f2,^) 

< \\Ur\\LHn^) \\DUr\\L3(f2^) II A^r 1 1 (fi,) 

+ c{\\DUr\\L!>{n^) + ||Ur||L3(r-2,)) \\DUr\\L^(n^) 1 1 A^Ur 1 1 (^-2^) 
1/2 II n„ l|l/2 



< c||Mr||L6(r-2,) \\DUr\\j^2i^i-2,) 1 1 -^^'^ 1 1 L« (fi,) \\^rUr\\L^(n,.) 

+ c| I £'Wr 11^(^(^-2^) \\DUr\\]^g^^j^^ ||^rUr||L2(r2,) 

+ c||Ur|li(^f-2^) ||Ur|li6^(f-2^)||£'Wr||L2(r2,) 1 1 ^r^r 1 1 (^2^ 



< c||w. 



1/2(12^) Il^'-|l_f/i(f2^) W^'^^Wn^in^) 1 1 -^"'-11/2(^-2^-) ||^rWr||L2(f2^) 
+ c(||£'Mr|li2^(^2,) + \\ArUr\\\(^f^i-2^^) 1 1 Du^ 1 1 ^(^(^^^^ \\ArUr\\ L^f2^) 



for almost every t £ [0,r]. We therefore conclude that 

^||V^.,|P + WArUrf < C(l + ||u.(0)||i2(M2))(l + ||Vzi,||i2(,2,))3. (40) 

Dropping || ArUr||?2f o ) and integrating, we obtain 



^ + \\^Mt)\\hif2^ 



< 



\Ur{0)\ 



2 



^1 - ct{i + h.(o)i|4,(^,^) (1 + wurmiiy 

Therefore for almost every t <T with 

^ ^ 2c(l + ||u,(0)||i,(^,^)(l + ||^,,(0)||2^,)2' (41) 
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we have 

Wurimmin) < ci\\urm\,\\urm\mAKm, koD- 

By Theorem [TJ and also since by ([55)1 . Af and J|a;p are uniformly 
bounded by ||Mr(0)||_f/^, we conclude that Ur is uniformly bounded in 
L°°{0,T;H^{M.'^)). With the above bound, we integrate (gO]) and obtain a 
uniform bound on \\ArUr\\L^(o^T;L^{Or))- This, by Lemma|4]and the uniform 
bound on \\ur{t)\\jfi(^j-2)j implies that 

Using ([Mt . we have 

/ \h\Ut<c sup Wurimmin^) [ hrWhin.^ < C. (42) 
Jo te[o,T] Jo 

Finally since iir = —ArUr — Pr/3(ur — hr, Ur), we have 

||w||l2(js.2) < c||ylrUr||L2(R2) +c\\{{Ur - K) ■ V)Mr|| L2(f2,) 
+ c|/lr| ||Ur(i)||ifi(f2,0 

which gives a uniform bound on \\ur\\L^(^Q^T:L^{M^))- O 

Lemma 6. Let i32\) hold and consider {ur,hr) to be a strong solution of 
with Ur and hr uniformly bounded in 

L'^{Q,T; H^ifir)) n L'^{0,T; H^{M.^)) n H^{0,T; L^iR"^)) and L'*(0,T;R2) 

respectively. Then there exists 

u e L^{0, T; n L°°{0, T; H^R^) n H^{0, T; L^{R^)) 

and h E L'*(0,r;K^) and a real positive decreasing sequence {r„}„gN con- 
verging to zero, such that, as n — > oo, converges strongly to u in 
L^((),T; H^{[2s)) for any s > 0, and hr^ converges strongly to h in 
i2(0,T;M2). 

Proof. By Lemma [3 u € L°°{0,T; H^{R^) n H^{0,T; L^{R^)). To show 
that u € L^(0, T; 7J^(R^)), we note that since Ur is uniformly bounded in 
L2(0, T; iJ^(J7r)), for any s > r we have 

\\ur\\L-^(a,T-H^{ns}) ^ ^ (43) 

with C independent of r and s. Therefore for any s > 0, there exists a 
sequence {ur„}„eN which converges weakly to u € L?{Q,T\ H'^{Qs)) as 
n — >■ oo. Hence, by uniqueness of weak limits u £ L^{Q,T] H^{Qs)) for 
any s > 0. Now we note that by Theorem 2.1.4 of Ziemer (1989) a function 
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w e L'^{Q) with C M'^, is in H^{f2) if and only if w is absolutely contin- 
uous on almost all lines in f2 parallel to the coordinate axes, and its partial 
derivatives belong to L'^{Q). Consider a real positive sequence {s„}„gN con- 
verging to zero as n — > oo. Since u E L^{{),T] H^{Qs^)) for any n, Du is 
absolutely continuous on all lines parallel to coordinates axes in ^nd 
the partial derivatives of all elements of Du are uniformly (in n) bounded in 
L'^{[2s^). Then, since n.^^iR^ \ I2s„ = {(0,0)}, and noting that the count- 
able intersection of sets of full measure in [0, T], has full measure in [0, T], it 
follows that Du e ^^(o, T; ^(R^)). This then, since u G £^(0, T; ^(R^)), 
implies that u € L^{0,T; H^{M.^)). 

We now show the strong convergence in L^{0, T; H^{f2s)) for any s > 0. 
Since Ur is uniformly bounded in L^{0, T; nL°°(0, T; Hr) n H^{0, T; V;), 
it is uniformly bounded in L^{0,T; 14)nL°°(0,r; H,)nH'^{0,T; V*) for any 
fixed s > r and therefore there exists a subsequence {ur„}neN (with r„ < rm 
if n > m) which converges strongly in L^{0, T; Hs) as n — >■ oo (see Robinson 
2001 or Temam 1977 for example). Therefore for any e > there exists 
some N such that for n > 

I!" - '«r„||L2(o_T;L2(R2)) < e. (44) 

For any fixed s > 0, let Es[u — be an extension of m — from to 
R^, such that for almost every t G [0,T] 

Es[u — Ur^] = u — almost everywhere in fis, (45) 

and 

\\Es[u-UrJ\\H''(R^) < Ck\\u- UrJ\Hi'(r2,)- (46) 

with fc > and = Cfc(s) (one can show that cq is independent of s but for 
our purpose here it makes no difference). By Ehrling's lemma for any 77 > 
there exists a c,, > such that 

\\Es[u - Ur,J||L2(o,T;ffi(R2)) 

< T]\\Es[u ~ Ur,J||L2(o,T:ff2(j!.2)) -|- Cjj\\Es[u - Wr„] || (0,T;L2 (r2)) . 

Let e > and take 77 — e/c2(s). For any n > m, we can write 

II" ~ '"r„||L2(o,T://i(r2,)) < 1 1 -E's " '"■r„] 1 1 (0,T;_f/i (R2 )) 

< 1]\\Es[u - Ur„]|lL2(o,T://2(R2)) + Crf\\Es[u - Mr„] || L2(0,T;L2(K2)) 

< e\\u - UrJ\L-^(Q^T-m{n,y) + Cq C^||w - UrJ\L^(Q^T■L■^i^2,))■ 

By for n big enough we have 

II II ^ 

CoCr, 

Also by ^ 

\\u - UrJ\L-2{0,T;H2{{2,)) <C 
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with C independent of n and s. Hence, for any e > there exists some Ni 
depending on e and s such that 

\W - Ur„ \\l^{0.T:HHO,)) < Ce 

for n > Ni. Therefore for a fixed s, ||u — Mr,, \\L^{o.T;H^{ns)) — >■ as n — > oo. 
It remains to show that /i^ — > A in L^{0, T; M^). For any n > m we have 

|(/ir,„ — /ir„) ■ n| 

= I ("r„, (fmll) - Ur„ (7'«n)) • n| 

= |(ur„(r™n) - Ur,^{r„in)) ■ n+ (ur,^{r„in) - ^^^(rnn)) • n| 

< sup |Mr„,(2^) -Ur„(2^)| +C||u,.„||if3/2(f2, 

<r II _ 11^/2 II _ 11^/2 

_ '^ll^r,,, ^^Ji II //I (^2,. ) ll^^m ^^Ti 1 1 //2 (■ J 

+ cIIm 11^/^ llii 11^/^ r 

^ ''ll"'-nllffi(r2,-„) Il"'""llff2(^^^) 'm 

with a < I. Since the above inequahty is true for any n, we conclude that 
for n > m > N , 

T 

\hr^ - ^r„ P dt -)■ as N —i' oo 

and obtain the required strong convergence. □ 

The above lemma shows the strong convergence of subsequences of Ur 
and hr- However, after showing that the limit u satisfies the Navier-Stokes 
equations in the whole plane and therefore is unique, one can show the above 
strong convergence result for any r rather than only a subsequence {r„}: 
we argue along the lines of the proof of Lemma 3.1 of Robinson (2004). 
We suppose that for some s > 0, Ur does not converge to u (strongly) in 
L'^{0,T; H^{f2s))- This means that there exists a subsequence {Mr„,} and a 
S > such that 

But , by Lemma [5] is uniformly bounded in 

L2(0, T; H^iQr)) n L°°(0, T; H^{m?)) n H^{Q, T; L^{R^)) 

and therefore by Lemma [6j has a subsequence that converges to u strongly 
in i^(0, T; H^{Qs)), contradicting the above inequality. One can argue sim- 
ilarly and conclude that hr ^ h in L^(0, T; R^). Therefore, from now on, we 
consider Ur (rather than Mr„) converging strongly to u in L^{0,T; H^{f2s)) 
for any s > 0, and hr (rather than /ir„) converging strongly to h in 
L2(0,T;M2). 

Corollary 1. Let Ur, hr, u and h be as in Lemma\^ Then 
l|u-Ur||L2(o,T;_f/i(^2r)) ^ 0' asr^Q. 
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Proof. Let < r < s. We can write 




< ||u - Ur||L2(o^T;Hl(f2,)) + C S \\u ~ 1 1 ^2 (o^T;H2 (f^,,)) 



using Holder inequality in the second line. Since 11" ~ 'itr||^2(Q 7^.^2(17 •)•) < C* 
with C independent of r, for any e > there exists a 5 (independent of r) 
such that s < (5 is small enough to ensure that 

cs ||u — u,-|||2(Q jn.^2(f-2^)) < e/2. 

Also, by Lemma |6] for r small enough 

The result therefore follows. □ 



5. The equations at the zero limit of the rigid body radius 

Now we can prove the main result of this paper: 

Theorem 2. Let i32\} hold and consider {ur,hr) to be a solution of @)- 
with the initial condition Ur{0) G H^{M?) satisfying Ur{0) = /ii-(O) for 
y e df2r and hr{0) G M^. Suppose also that there exists uq £ H^{M.'^) such 
that Ur{0) Uq in L^(M^) as r — > and \\ur{0)\\H'^{f2r) — c||uo||_ffi(R2) for 
any r > 0. Then 

i) u, the limit of Ur at the zero limit of r satisfies the Navier-Stokes equa- 
tions in the whole plane with initial condition uq, 

ii) the trajectory of the centre of the rigid disc converges to a fluid particle 
trajectory as r 0. 

Proof, i) We first show the result for t G [0,T], with T as in (gT]), and 
then use the properties of the solution of the Navier-Stokes equations in the 
whole of to get the result for any < t < 00 . In what follows 

H ^{u£ [L^{M^)f : V • u = 0} 



and 



V = {ue [i?i(R2)]2 : V • w = 0}. 



Limiting motion of a fluid-rigid body system 



21 



V* is the dual of V, and 



{u,v)n, = / u-vdy. 



1) {ArUr - Au) -vdydt-^-O as r for any v e L^(0, T; i/): 

We show the weak-* convergence of ArUr to Au in -L^(0, T, V*) and then 

use the fact that \\ur\\ l^(o,T;H^ (0^)) ^^'^ \ l'^ {o,t-h^ {Qr-)) bounded inde- 
pendent of r, to conclude the required convergence for and v G L^(0, T; H). 
Let V e i^(0, T; t/) and Vr = PrV and write 



{ArUr — Au,v) dt — / (yl^Ur, Wr) -I- / {l'Au,v)dt 



{ArUr + l^Au,v) dt + / (yl^Wr, Wr — dt. 

(47) 

We have 

{ArUr,Vr - v) dt < C \\Ur\\ (Q^T-H'' in^))\\'"r " v\\ (^q^T;L^ (n,)) 







{vr — v) dy dt 



+ / 4 / D{ur)ndr- [ 
Jo Jonr J B, 

+ [ [ ■ D{ur)ndr [ {vr - v) ■ y^ dy dt. 

Jo ^ Jdi2r J 



By characterization of G2 in Lemma [TJ the last term is zero. For the second 
term in the right-hand side of the above inequality, by Lemma [5] and [3] we 
can write 



D{ur)ndy ■ j {vr — v)dy 



^ Id ^'<^^^'^'^y ' j {vr-v)dy^ 

< -J \\D{ur)\\L-ion^) (27rr)i/2 \\v - Vr\\L-(B^) (^r2)i/2 

< Cr"~^/^ \\Ur\\H^f2^) l|l'llifi(R2)- 

for any a < 1, and therefore 

/ {ArUr,Vr — v) dt ^ aS T — > 0. 

Jo 
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For the integrand in the first term of the last line of (j47p we can write 



{ArUr + vAu, v)=-2v I {Diur) - D{u)) : D{v) dy 

J sir 

+ 2v j {D{ur) - D{u))v ndr + Au-vdy 

-T7 [ / D{ur)ndr 

A/ Jb,. Jdn,. 

— -y^ I V ■ dy j ■ D{ur)ndr. 

J Jb,. JdB,. 

Since with < a < 1, 

|- / vdy [ \D{ur)dr 
Jb,. Jan,. 

< -\\vU^B,.){^r^f'''\\Dur\\ 



—r j V ■ y^ dy n'^ ■ D{ur)ndr 

^ ^Ml' (^')'^' \\DUr\\mdB,.)V2^ 
< Cr"/2+l/2 ||w||ffl(R2) \\Ur\\H2(n,.), 

and 

/ Au-vdy < \\u\\H2(m2)\\v\\L4(^B,.)\Br\^^'^ 

JB,. 

we have 

/ {ArUr^ vAu,v)dt < C || U - || (o,T:Hi (R^)) 1 1 1' 1 1 (0,T:if i (R^)) 
JO 

+ Cr^WUr - m||l2(o^T;H2(r2)) \\v\\l2(o^t-H^{R2)) 

+ Cr^/^ ||Ur||L2(0,T;//2(R2)) |jtl||i2(o .r.^l(R2)) 

which converges to zero as r tends to zero. We conclude that Aj-u,. — Au ^ 
in L2(o,r;F*). 

Now to prove the required convergence for v E L^{0, T; H), we note that 
since L^{0,T;V) is dense in L^{0,T;H), for any w e L^{0,T;H), there 
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exists a sequence {vm}m=i -^^(0) ^)) such that Vm ^ w in 1/^(0, T; H) 



as m — > 00. We can write 



Jo 

/ {ArUr - Au, W - Vm)n^ dt + {ArUr - Au, Vm)Qr 

Jo Jo 



< \\ArUr - AM||L2(o,T;L2(r2r))||w " «m 1 1 L2 (o,T;L2 



+ 



/ {ArUr - Au,Vm)nr 

Jo 



dt 



For any e > 0, there exists M big enough such that \\w — Vm\\ < e/{2C) 

for m > M and with C > \\ArUr — ^w|lL2(o.T:_tf(r2r))- Also by the weak-* 
convergence we just showed, there exists some r = r{M) such that 

/ {ArUr - Au, Vm)nr dt < e/2, for r < r(M). 
Jo 

Therefore for any e > 0, 

/ {ArUr — Au, w)q^ dt 

Jo 



< e 



when r is small enough. 

2) fr\l3{ur - hr,Ur)\ V[{{u - h) ■ V)u] in L^{0,T;H): 

We note that P[((u - h) ■ V)u] is bounded in L^(0,T; H), since 

< \\{U- V)M|||2(o,T;iJ) + V)u|||2(o,T;Jf) 

< l|w||L~(0,T;ffi(R2)) ll"llL2(o,T;i?2(M2)) 

+ I|/i||l''(0,T;R2) ||w||i,4(o,T;H1(R2))- (48) 

We first show the weak-* convergence for ^^/^(O, T; F*). Let u e i^(0, T; if^) 

and Vr — PrV as before and write 

/ {Fr[P{Ur - hr,Ur)] - P[((w - fl) ■ V)u],v) dt 

Jo 



{[{Ur — hr) ■ VjWr, Vr 

Jo 

{[{Ur - hr) ■ S/]Ur - [{u - h') ■ S/]u,v) ^ dt 
+ / {[{Ur - hr) ■'^\Ur,Vr - v) ^ dt + {[{u ~ h) ■ V]u, v) dt. 

Jo ''Jo 
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We have 

(■T 



I (((Wr — hr) ■ V)Ur,Vr " ^ dt 

Jo 

< (^||Wr||l,~(0,T;/J'i(f2^))||Wr||z,2(o,T;H-2(r2^)) 

+ ||^||z,2(0,T;R2)||Wr||L°°(0,T;ffi(r2,))) ll^" " lU^ (0,T;L2 (r?, )) , 

/ {{hr ■ V)Ur - {h ■ V)m, v) „ dt 
Jo 







{{{hr - h) ■ W)ur, v)^^ + {{h ■ V)K - w), v)^^ dt 

<\\hr - /l||z,2(o,T;R2) ||Wr||L~(0,T;H-i(f2^)) ||i^||l2(0,T;L2(k2)) 
+ ||^||z-4(0,T;IR2) \\Ur - w||L2(0,T;/J'i(f2^)) ||^^|U4(0,T;L2(k2)) , 



with W = Ur — u 

i-T 



I {{Ur ■V)Ur — {u-W)u,v) 

Jo 

= [ {{wy)w),v)^ +{{wy)u),v)^ +{{u-V)w),v)^ dt 
Jo 

< C (||w||i:,oc(o,T;ffi(r2,)) + ||w||z,~(0,T;ffi(^2.))) 

X \\w\\L^{0,T;m{nr)) \\v\\L^{0,T;m(nr))' 

and 

/ {[{u-h)-V]u,v) dt 
Jo 

< J C (^||u|jHi(K2) ||u||^2(K2) + |/l| ||w||i4(B^) jB^I^/" 

< (||/l||z,2(0,T;R2) + ||w||L2(0,T;ff2(M2))) 

X ll'"IU~(0,T;Hi(M2)) ||l^||L2(0,T;fi'i(M2)) 

We therefore conclude the required weak-* convergence. The weak conver- 
gence in L^{0, T; H) then follows in a similar way to the previous case, since 
both PJ/3( )] and P[((u — h) ■ V)u] are bounded in L'^(0, T; i?) and 

i^(0, T; is dense in ^^(o, T; il). 

3) /(f (dK - w)/dt, t-)^ dt^O as r ^ for any w e L^{Q, T; H): 
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As before, we first show the convergence for more regular v e i^(0, T; V). 
For s > r we can write 



^ / du 
dT ~ dt ^ 

^ / dur du 
dT ^ di 



For V € L2(o^ we have 



dt 








) dt + 


n 







dur du 

"dT ~ dt ' 



dt. 



dur du 

"dT ~ dt ' 



dt 



< 



f 


dur du 


Jo 


~dr " dt 



\^2r\n,\'^^\\v\\LH^2Ans} 



< CS^/^||u||i2(o^T;V') 

since \\ur\\L'^(o^T:L^(nr)) uniformly bounded by Lenima[5]and |1u||l2(q 7-.2;,2(r2)) 
is bounded by Lemma [6) Therefore, for any e there exists S > such that 
for any r and s less than 6 



/ dur du 
V"dt~~dt'^ 



dt < e/2. 



By Lemmaini u^ ^ u as ?- ^ in L^(0, T; L'^{f]s)) for any fixed s and hence 
there exists r* > such that 



/ dur du , , 



for r < r* . Let r — min{S/2,r*}. Then 
rT 



dur du 

-r^--r,v] dt<e 
dt' 



and we get the required convergence for v G L^{0, T; V). 

The converegence for v E L^{0,T; H) follows in a similar way to the 
parts (1) and (2), since \\ur\\L'2{o,T;L^inr-)} and ||u||i2(o,T;L2(R2)) are both 
bounded by a constant independent of r. 

The above convergence results and (fT9|) imply that u(t) satisfies 

'■'^ ^du 



/O ^f2r 



dt 



Ai 



- h) ■ V)u] • w dy dt -> as r -> 0. 
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for any v G L'^{0,T;H). Since u,Au e L'^{0,T;H) by Lemma E] and 
P[((-u - h) ■ V)m] e £2(0^ T; H)hy in the zero limit of r, we have 

+ Au + ¥[{{u~ h) ■V)u] ^0 as an equahty in L2(0,r;i7). (49) 



We need to show the limiting velocity field u has initial condition uq. 
Let (j) e (7^(0, T;i?) with (/)(T) = 0. Taking the inner product of the above 
equation and (jl9p with (f>{t) we obtain 

/ (u, ^)n^ + (Au + P[((u - h) ■ V)u],q^)t2,. dt = (u(0), (/)(0))fi. 

and 

[ {Ur, + {ArUr + P.[/3(w - h, u)],C^)a^ dt = {Ur{0), 0(O))f2, 

Jo 

respectively. Comparing these two at the zero limit of r implies that u{0) = 

Uq. 

Since h e H^{0,T;R'^) C L°°{0,T;M.'^), we can perform the following 
change of coordinates in (|49|) 

x{t)=y{t) + h{t). 

Letting 

U{x,t) ^u{x + hit),t) 

we can write (|49)) as 

^ +AC/ + P[(C/- V)[/] =0 as an equality in £2(0, T;F). (50) 



Now, we extend the result to all times: 

Taking the inner product of the above equation with AU, using appro- 
priate estimates for the nonlinear term and then integrating, we obtain 

||u||hi(r2) ce=ll"«ll^^ ||"o||hi(r2). (51) 

Let ||Mr(0)||Hi(j7^) = C\. By Lemma[5l we have ||ur(t)||^fi(f2r) < 2Ci for 

*<^*^ 2c(l + ||«,(0)|||,J(l + C2)2 < ^ 
(noting that ||wr(Ollffr — ll"i-(0)||_ff,, for any t > 0). Since 



T 



\\ur — ^0 as r 0, 
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there exists r* small enough such that for some t > T ~ T* /2 

\\ur>{t)-u{t)\\H.^^n^,)<Ci/2 

and therefore |lur* (i)||ffi(r2^.) < Ci and by above for t < T + T*/2, 
\\ur{t)\\ (n^) is uniformly bounded. Hence, by the first part of the proof 
u satisfies ^ for t € [0,T + T*/2]. The same estimate ([5T|) holds for 
0<t<T + T*/2as well. Therefore the same argument as above gives the 
result of the theorem for < t < T + T* and continuing in a similar way 
for < t < T + nT*/2 for any n e N. 

ii) To show the convergence of the trajectory of the centre of the tracer, 
we fix a point y* on df2r and write 

|u((0,0),t)-/i(i)| < HiO,0),t)~u{y*,t)\ 

+ \uiy*,t)-hr{t)\ + \hrit)~hit)\. 

We have 

\u{{0,0),t) - uiy*,t)\ < cr"/^ h||ff3/.(R.) 

for a < 1 by Lemma 7.26 of Renardy & Rogers (2004). To find an appropri- 
ate bound on the second term in the right-hand side of the above inequality, 
we note that for any n = y/\y\ with y e df2r we can write 

\n-iuiy*,t)~hrm 

= |n- iu{y*,t) - Ur{y,t))\ 

< \u{y*,t) - u{y, t)\ + \u{y, t) - Ur{y, t)\ 

< Cr"/^ ||u||ff3/2(R2) +c\\u-Ur\\]^^|^|-^^^ 

Since the above inequality is true for all n as above, we conclude that 
\u{y*,t)-hr{t)\ <cr"/2 ||u||^3/,^^,j 

, II „ ||l/2 II _ ||l/2 

Therefore we have 

r\uiio,o),t)-hit)\^dt 

Jo 

< Cr" ||u||L~(0,T;ffl(B2)) |1"|1l2(0,T;H2(r2)) 

+ c\\u- Ur\\L-2{Q,T:H2{{2^)) \\u - 1 1 (o,T;Hi (^"^r)) + ll^'' ^ ^11 L2(0,T;B2) ■ 

The right-hand side converges to zero as r — > by Lemma |6l Hence 

[ \h{t)-U{h{t),t)\dt = 0, (52) 
Jo 

implying that h{t) is the trajectory of the fluid particle initially at (0,0). 
Since 

\h{t) - hr{t)\ < \\h- hr\\miO,T:R^), 

we conclude that hr{t) h{t) as r -> by Lemma (O • □ 
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We note that there exist uq and Ur(0) satisfying the assumption of The- 
orem [2l For instance we can consider uq G If^(R^) such that uq = h for 
y G Bn with R> and some /i G M^, and let ^,.(0) = mq for any r < R. 

6. Conclusion 

We studied the limiting motion of a system of a rigid disc moving with 
a fluid flow in as the radius of the disc goes to zero. We showed that 
if the disc is not allowed to rotate, the trajectory of the centre of the disc 
converges to a fluid particle trajectory. 

Two related problems are the motion of the fluid-rigid body in a two- 
dimensional bounded domain, and the case of more than one rigid disc. 
In the case of a bounded domain, the existence of a global strong solution 
of a fluid-rigid body system is shown by Takahashi (2003) assuming that 
the rigid body does not touch the boundary. With the same assumption, it 
seems that a similar approach as the one presented here (but more technical 
since the change of coordinates would be more complicated) can give the 
same convergence result. When there are several rigid bodies moving with 
the fluid, which is the case in real experiments, it is known that at least one 
weak solution exists when the radius of rigid bodies are fixed (Desjardins & 
Esteban 1999; Feireisl 2002; Grandmont & Maday 2000; San Martin et al 
2002), but to our knowledge, the uniqueness of the solution is not known. 

The analysis here goes some way towards justifying the use of tracer 
particles for finding Lagrangian paths of fluid flow, but we have unfortu- 
nately had to restrict to the case in which the particles are not allowed to 
rotate. Obtaining appropriate uniform bounds in the case that the parti- 
cle can rotate seems more challenging (see Remark [T]), and presents a very 
interesting open problem. 
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